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In this work, we introduce long range version of the extended Hubbard model. 
The system is defined on a non-uniform lattice. We show that the system is inte- 
grable. The ground state, the ground state energy, and the energy spectrum are also 
found for the system. Another long range version of the extended Hubbard model is 
also introduced on a uniform lattice, and this system is proven to be integrable. 
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In recent years, there have been considerable interests in low dimensional electronic 
systems of low dimensions. Systems of this type may exhibit interesting novel physics, due 
to low dimensionality and strong correlation. Anderson has suggested that the 2D one band 
Hubbard model should explain the basic physics of the high temperature superconductivity 
It is suggested that the normal state of the new cuperate oxide superconductors may 
share the feature (Luttinger-liquid-like) of ID interacting electron gas This physical 

motivation has initiated lots of recent activities in one dimensional electronic models. 

In one dimension, the Hubbard model was exactly solved by Lieb and Wu With 
the 5*0(4) symmetry due to Yang and Zhang the completeness of the eigenstates 

in this structure has been discussed 0. More recently, an extended Hubbard model was 
introduced by several authors in general d dimensions [^]. By using r^-paring mechanism, they 
construct eigenstates of the system, with off- diagonal-long-range order 0. Particularly, its 
one dimensional version was solved exactly with Bethe-ansatz . In this work, we introduce 
the long range version of the extended Hubbard model in one dimension. We show that the 
system is completely integrable, with explicit construction of infinite constants of motion. 
The ground state, the energy spectrum are found for the system. 

Let us consider a one dimension lattice of sites L. The positions of the sites are given 
by the roots of the Hermite polynomial Hl{x). The roots of the Hermite polynomial 
'"i ; ''"2 5 ■ ■ ■ 5 ''"-L are all real and distinct, and this non-uniform lattice is thus well-defined. 
Consider the following Hamiltonian: 

H = -l E J^^^^v (1) 
where Jij = l/(rj — r^)^, and the permutation operator Hj^ is given by 

Uij = c]|Qt(l - riii - riji) + c\^Cj^{l - rin - rij^) 
+cj^Q|(l - - rij^) + 4^9^(1 - nil - rij^) 

- - 1) + clAl^jlCn + c^CiTsVii 
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+ - l/2)(n,i - 1/2) + (n,^ - l/2)(n,; - 1/2), (2) 

where Ci„ and c|g. are electron annihilation and creation operators at site i with spin cr =1, |. 
The electron number operators are given by rii^ = c\-^c^,nii = cl^Cii,ni = rii^ + nj|. If one 
defines the Hamiltonian on a uniform lattice and the coupling parameter Jij = 6i^\i-j\, the 
system becomes the extended Hubbard model studied by Ebler, Korepin and Schoutens 0. 

For the long range extended Hubbard model defined on this non-uniform lattice, we'll 
present its ground state, the energy spectrum and the proof of integrability in the following. 
At each site i, one can use the standard slave boson representation: 



itxO 


= plh 


iT><^ 


= Pifia 


cr >< r 


f icr fir 


a >< 0| 


= flh 


>< a\ 


= blU, 



where | iT>= cJ|C-|^|0 >, \a >= cl^\0 >, at site i. For these relations to hold, the pair 
operators p and the hole operators b are bosonic, while the operators / are fermionic: 

[bi,b]] = 6ij 
[Pi,p]] = Sij 

[fia,fl]+ = Sij6^r, (4) 

with the constrain that at any site i there is always one particle, i.e., J2a=i,i fLfi^r + + 
pjpi = 1. With the slave boson representation, a state vector can be written as below: 

A M Q 

■I[pLUfLIibl\0>, (5) 

«=1 1=1 i=l 
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where the wavefunction is symmetric when exchanging {x} or {z} respectively, while 
anti-symmetric when exchanging ?/j(Tj and yj(Jj. In the following, we use {qi,q2, ■ ■ ■ Al) = 

{Xi, - ■ ■ ,XA\yi, - ■ ■ ,yM\zi, - ■ ■ ,Zq). 

With the Hamiltonian H given by Eq. (|I]), the number of the electron pairs, the number 
of holes, the number of the sites single-occupied by up-spin electrons, the number of the 
sites single-occupied by the down-spin electrons, are all conserved quantities. Let us denote 
them by the following notations: 



i=l 
L 



i=l 
L 



i=l 

Q = L-M^-M^-A 

M = M^ + M^. (6) 

In the following, we work in the Hilbert space of fixed A, M^, M| and Q. 

Using the amplitude (p defined by the Eq. (H), the eigenenergy equation can thus reduce 

to 

-I E 7-^M,,0({g}) = (7) 

where the operator Mij permutes the coordinates and qj. The Hamiltonian H = 
~| J2i<i^j<L{qi ~ (Ij)^'^^ij commutes with the infinite number of simultaneous constants of 
motion, i.e., [H, In] = 0, [/„, /„] = 0, with n,m = 1,2, oo, where /„ = J2i=i = '^l'^*; 

with a] = Efc(^j)=i i{qj - qk^^Mjk + iqj, ai = {a^Y , with j = 1, 2, ■ ■ ■ , L This shows 

that this long range version of extended Hubbard model is also integrable. When there 
are no electron pairs, this system reduces to the long range t-J model studied before RITU 



When there are no electron pairs and no holes, the system becomes a previous spin chain 
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We first construct the ground state wavefunction in the general subspace of fixed non-zero 
A, M|, M| and Q. Let us look at the following Jastrow product wavefunction: 

(j){Xi,X2, ■ ■ ■ , XaIViCTi, y2<y2, ■ ■ ■ , VmCTmIzi, ^2, " " " , Zq) 

= n (l/i-l/j)^'''''^e^5"^'"(""-"^). (8) 

l<i<j<M 

This wavefunction is anticipated to be the ground state of the system in the subspace of 
fixed A, M|, and Q. With previous experience |10|, one may show that this wavefunction 



is indeed an eigenstate of the Hamiltonian, with the following eigenenergy 

Eo = -\l{L - 1) + Im^Mt - 1) + \m^{M^ - 1). (9) 

For this Hamiltonian, one may show that a] and aj are raising or lowering operators 
defined on this lattice, [a], H\ = —a], [aj, H] = aj. One may simply prove following identities: 

0^0 = 0; 2 = 1,2, ■■■,M 

0^0 = 0; i = A + M+l,A + M + 2,---,L. (10) 
Therefore, one can see that 

M 

($:a,)0 = O, 

i=l 

( E a,)0 = O. (11) 



i=A+M+l 



Furthermore, one can show that 



M 

(Ea.+Aan0 = O. (12) 

i=l 

The above three identities show the impossibilities of constructing non-zero eigenstates of 
lower energy than Eq with the annihilation operators. One may regard this as a partial 
support for our idea that (j) is the ground state in the subspace of fixed A, M|, M^, Q. For a 
more complete confirmation, further work is needed in this case. 

The full energy spectra of the system consists of energy levels equal-spaced. 



E{s) ^Eo + s, (13) 

where s = 0, 1, 2, • • •. For a lattice of finite size L, there is an upper bound on tfie value of s, 
as the Hilbert space of the system is finite. There are several ways of constructing excited 
states. The first way is to excite those electron pairs: 

|ni, n2, • • • , >= E Ui^'lT"' 1'^ (14) 
p i=i 

where ni,n2, ■ ■ ■ ,nA are integers or zero, the summation P is over all possible permuta- 
tions. Under the operation of the creation operators, any states constructed this way, if not 
vanishing, will be the eigenstates of the Hamiltonian with eigenenergy 

Ei^ Eo + {ni+n2^ Ktia). (15) 

The second way to create excitations is to excite the holes of the system: 

p 

with the eigenenergies given by 

E2 ^ Eo + {mi + 1712 + ■ ■ ■ + mg), (17) 

if these state vectors are not zero. Finally, one may create excitations by exciting the 
electrons on the sites single-occupied. These excited states may be written as 

103 >= E('5l(^l)^)""K'52(i^2))^)""^ ■ • ■ (^M(i^M)^)""-|0 >, (18) 
P 

where Sil^iY = a|^_,_^(T|'% with i/i = 0,±, z. These states have eigenenergies given by 

Es = Eo + {ni + n2 + --- + nM)- (19) 

Clearly, excited energy levels are highly degenerate. However, we are still unable to de- 
velop a systematic rule to characterize the pattern of the degeneracy and to explain it with 
symmetries of the system (presumably, Yangian symmetry) . 
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In the limiting case where M = + Mi = 0, one has two types of bosons on the chain, 
the local electron pairs and the holes. The wavefunction becomes 

X2, ■ ■ ■ , xa\zi, Z2,---, Zq) = 1, (20) 

which is obviously the ground state of the Hamiltonian H, with the ground state energy 
Eq = —L{L — l)/4. With the raising operators, one may construct the first excited state 
quite easily: 

(l)'{xi,X2, ■ ■ ■ , Xa\zi, Z2, - ■ ■ ,Zq) = {Xi + X2 ^ h Xa)- (21) 

The wavefunctions of higher excited states may be obtained by acting this wavefunction 
with raising operators again. 

Besides this extended Hubbard model of long range type defined on a non-uniform lattice, 
we would like to introduce another model of extended Hubbard model of long range type. 
Consider a uniform lattice of sites equal-spaced. The size of the lattice is L. The Hamiltonian 
of the system is given by 

H = -l E (22) 

where the function d{i — j) = ^sm{n{i — j)/L), and the permutation operator H^- is the 
same as Eq. (Q). With our previous experience, one can simply write the first quantized 
Hamiltonian as follows: 

H = -l E l/d\q^-qJ)M,,. (23) 
It is very easy to prove that this extended Hubbard model is also integrable, once the 



Hamiltonian is written this way by using the relation of Fowler and Minahan [|1| 

L L 
Ln = E = E ~ Zj)Mij, 

[L„, Lm] = 0, [L„, H] = 0; n, m = 1, 2, • • • , (24) 

where Zj = e^'^*''^/^, j = 1, 2, ■ ■ ■ , L. This thus provides a proof for the integr ability of the 
extended Hubbard model. In the limiting case of no electron pairs, this system reduces 



to the long range t-J model [1T3H15|,^. If there are no holes and no electron pairs, the 



Hamiltonian reduces to the Haldane-Shastry spin chain |TT],[T2|. In the general case, to write 
the conserved quantities of the system in terms of second quantization, one may simply use 
the permutation symmetries of the wavefunction. In this paper, we shall not present the 
wavefunctions and the eigenenergy spectrum of this system. More results on this system 
will be given elsewhere. 

In summary, we have introduced two integrable models of strongly correlated electrons, 
one defined on a non-uniform lattice, the other on a uniform lattice. The systems belong to 
the Jastrow-integrable-type. Particularly, the system on the non-uniform lattice has equal- 
spacing energy spectrum. Because of the degrees of freedom of the paired electrons, the 
energy levels will have much larger degeneracy than the long range t-J model. One of the 
most interesting things is to determine the underlying symmetries explicitly that account 
for the degeneracy of each energy level. 

I wish to thank V. E. Korepin for interesting conversations. Discussions with Mo-lin Ge 
and C. Gruber are acknowledged gratefully. This work was supported by the Swiss National 
Science Foundation. 
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